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Abstract 

We present a formalism to evaluate QCD diagrams with a single vir- 
tual gluon using a running coupling constant at the vertices. This method, 
which corresponds to an all-order resummation of certain terms in a per- 
turbative series, provides a description of the momentum flow through the 
gluon propagator. It can be viewed as a generalization of the scale-setting 
prescription of Brodsky, Lepage and Mackenzie to all orders in perturba- 
tion theory. In particular, the approach can be used to investigate why 
in some cases the "typical" momenta in a loop diagram are different from 
the "natural" scale of the process. It offers an intuitive understanding of 
the appearance of infrared renormalons in perturbation theory and their 
connection to the rate of convergence of a perturbative series. Moreover, it 
allows one to separate short- and long-distance contributions by introduc- 
ing a hard factorization scale. Several applications to one- and two-scale 
problems are discussed in detail. 
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1 Introduction 



Perturbative expansions in field theory are asymptotic, making it necessary from a 
conceptual point of view (if not from a practical one) to truncate any perturbative 
series at some finite order in the coupling constant. This truncation introduces 
renormalization-scale and scheme dependences. In QCD, this poses a problem 
because of the strong scale dependence of the running coupling constant ^^(/i^) in 
the low-momentum region. In many cases this puts limitations on the precision 
with which experimental data can be described in terms of the fundamental 
theory of strong interactions. 

Several prescriptions on how to fix the renormalization scale and scheme in 



a truncated perturbative series have been proposed IH-lT^l- They all rely, in 
some way or another, on a guess about the size of uncalculated higher-order 
contributions. This guess can be based on criteria such as the apparent rate of 
convergence of a series, the size of the coefficient of the last term in the truncated 
series, or the sensitivity to changes of the renormalization scale and scheme. It 
can also rely on physical criteria such as the role of various mass scales in a given 
problem. An example for such a physical scheme is provided by the scale-setting 
prescription of Brodsky, Lepage and Mackenzie (BLM), which amounts to absorb 
certain vacuum polarization effects appearing at two-loop order into the one-loop 
running coupling constant ||^. Extensions of this scheme beyond the one-loop 
order have been considered in Refs. |]I0|-|]T2| . 



A large higher-order coefficient in a perturbative series can arise from an 
anomalously large contribution of a particular set of higher-order diagrams. For 
instance, one can imagine a series with a large two-loop coefficient, but small one- 
loop and higher-order coefficients. Such "genuine" higher-order effects are hard 
to anticipate without a detailed calculation. On the other hand, it is possible 
that large higher-order corrections result from an inappropriate choice of the 
renormalization scale or scheme. In particular, a change of the scale in the running 
coupling constant at one-loop order leads to a change proportional to Pq at 
two-loop order, where /3o = 11 — ^Uf is the first coefficient of the /^-function, 
and Uf denotes the number of light quark flavours. Since /3o is large in QCD, 
such effects can be numerically signiflcant. An appropriate choice of scale should 
try to minimize this type of higher-order terms. Surprisingly, however, one flnds 
that even in one-scale problems using the "natural" scale in the running coupling 
constant at one-loop order often leaves large corrections of order /?o«s- Let us 
consider some examples related to heavy quark systems: (i) the relation between 
the pole mass of the bottom quark and the mass rrib{mfj) renormalized in the 
modifled minimal subtraction (MS) scheme; (ii) the ratio of the decay constants 
of pseudoscalar and vector mesons, B and B* , in the so-called static limit; (iii) 
the parton model prediction for the semi-leptonic decay rate V{h ueve)- At 
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two- loop order, the corresponding perturbative series are [0]-|T5| 
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where the constant k is yet unknown. In all cases the running coupling constant is 
renormalized in the MS scheme. For n/ = 4, the parts of the two- loop coefficients 
proportional to Pq are 13.02, 4.40, and 26.08, respectively. At least in the ffist 
two examples these terms provide the dominant part of the two-loop corrections. 
Moreover, in all three cases the two-loop corrections are quite large, casting doubt 
on the convergence of the perturbative series. An important example related to 
light quarks is the so-called Z)-function, which is proportional to the derivative of 
the correlator of two vector currents containing massless quark fields. At two-loop 
order, one finds 
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where denotes the euclidean momentum transfer. Again, the term propor- 
tional to dominates the two- loop coefficient. However, the absolute size of the 
two-loop correction is smaller than in the cases considered above. 

The appearance of large corrections proportional to /3o can be interpreted 
as an inappropriate choice of scale in the one- loop running coupling constant. 
BLM have argued that one should absorb these terms, which arise from self- 
energy corrections to the gluon propagator, into the one-loop running coupling 
constant. This prescription defines the so-called BLM scale /xblm- For the above 
examples, one finds (in the MS scheme): hblm = 0.10 mf,, 0.21 mf,, 0.07 mj,, and 
0.71\/Q^. Clearly, if the BLM scale is to be interpreted as a "typical" scale of 
virtual momenta in the corresponding Feynman diagrams, the question arises 
why it is often much lower than the "natural" mass scale in the problem at hand. 

To analyse this issue, we propose a generalization of the BLM prescription. 
Consider the perturbative calculation of a physical (i.e. renormalization-scheme 
invariant and infrared finite) quantity S{M^), which depends upon a single large 
mass scale M. We restrict our discussion to Green functions without external 
gluons. A generalization of the method to physical cross sections and inclusive 
decay rates, which receive both virtual and real gluon corrections, will be pre- 
sented elsewhere |ff9l . The aim of the BLM prescription is to "guess" the average 



virtuality of the gluon in a loop diagram and to use it as the scale in the run- 
ning coupling constant. Clearly, a better way to proceed would be to perform 
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the calculation with a running coupling constant as{—k^) at the vertices, where 
k is the momentum flowing through the virtual gluon line. The result of such 
an improved calculation, which we denote by Sres{M'^) since it corresponds to a 
partial resummation of the perturbative series (see below), may be written as 



where f{k, M, . . .) is the integrand of the Feynman diagrams. The function w{t) 
describes the distribution of virtualities in the loop calculation. The integral 
over this distribution function corresponds to an average of the running coupling 
constant over the loop momenta in one-loop diagrams and, in a way, provides 
the optimal improvement that can be achieved without a complete higher-order 
calculation. 

The fact that the integration in (j^) extends to t ^ indicates the appearance 
of non-perturbative effects. It makes explicit that any perturbative calculation 
receives long-distance contributions from the integration over low momenta in 
Feynman diagrams. We shall discuss the significance of these contributions at 
length below. Here we just note that our approach goes beyond perturbation 
theory since it is equivalent to the resummation of an infinite number of terms 
in the perturbative series for the quantity S{M^). Let us write this series in the 
general form 



The expansion coefficients Sn can be written as a power series in the number 
of light quark fiavours or, equivalently, as a series in powers of Pq. We restrict 
ourselves to cases where the calculation of the one-loop coefficient Si does not 
involve the non-abelian gluon self couplings, and where there are no external 
gluons. Then the first dependence on nf comes at two- loop order, and in general 
the coefficients can be written in the form 



Sn = CnPr' + dnPr' + ■ ■ ■ + K = {' fn^)"" + 0(71]-^) . (5) 



The coefficient c„ in front of the highest power of f3o is proportional to the coef- 
ficient in front of the highest power of nj. This implies that c„ can be calculated 
from "quasi one-loop" diagrams, in which a gluon line is dressed by [n — 1) light- 
quark loops. It is the partial series built up by the terms proportional to c„ which 
is resummed in our approach. To see this, we use the one-loop renormalization- 
group equation for the running coupling constant to relate ^^(tM^) in to 
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Under the above assumptions it is easy to see that the running of the couphng 
constant is the only source of terms of order P^'^a^ in the series @. We thus 
conclude that 

The BLM prescription is to absorb the two-loop term C2 /3o «s i^^o ^ redefinition 
of the scale used in the coupling constant at one-loop order. It is equivalent to 
setting wblm(^) = Ci5{t — e"'^^/'^^), thus choosing an average virtuality /iBLM = 
exp(— C2/C1) for the gluon. Already at this point it is clear that this can only 
be a good approximation if the distribution function w{t) is narrow. If it is wide, 
it is better to perform the integral in which resums all terms of the form 

Cn /3q tt". 

Clearly, the resummation in (|^) does not replace a full higher-order calcula- 
tion. For instance, "genuine" two-loop corrections not related to the running of 
the coupling constant are not taken into account. Nevertheless, our approach can 
be considered as an optimal improvement of one-loop calculations, which takes 
into account the full information contained in one-loop diagrams combined with 
the running of the coupling constant. As such, we believe that the construction 
of the distribution function w{t) is an interesting new concept, which provides 
information that goes beyond what is contained in a low-order perturbative cal- 
culation. 

In Sect. we generalize the above discussion in such a way that the scale and 
scheme independence of the procedure become apparent. In Sect. ^ we discuss 
the asymptotic behaviour of the distribution function and its relation to the 
ultraviolet and infrared properties of the perturbative series. In particular, we 
trace the appearance of infrared renormalons and discuss how to separate short- 
and long-distance contributions by introducing a hard factorization scale. In 
Sects. I and ^ we calculate the distribution function for several one- and two- 
scale problems, among them most of the quantities considered in (|1|) and (Q). 
We present a detailed numerical analysis, in which we compare the results for 
the resummed series to low-order calculations and investigate the relative size 
of short- and long-distance contributions. Section ^ contains a summary and 
conclusions. 

2 Construction of the distribution function 

We will now repeat the above argument in a slightly more general form, which 
allows us to keep track of scale and scheme dependence. Consider some dimen- 
sionless, infrared-safe and renormalization-scheme invariant quantity S{M'^,z), 
which depends upon some large mass scale M and, in the most general case, on 
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a set of dimensionless parameters z. For instance, in a two-scale problem we may 
choose = m\m2 and z = 7712/ mi. Let us investigate the perturbative series of 
S{M'^,z) in powers of the coupling constant af{fi'^) renormalized at some scale 
H and in some renormalization scheme R: 



5pert(M^.) = f:(^p„^ 



M, z) . (8) 

The scale and scheme dependence in (||) cancels between the coefficient functions 
and the coupling constant. As in (^), the coefficients can be expanded in powers 
of po, 

S^ifi, M, z) = c^(/i, M, z) (3^-' + ..., (9) 
and we define a function wnit, yU, M, z) such that 



n-l 




S,UM',z) = Jdtwn{t,f,,M,z)^^Y^ = j:l^^) 

Q n=l \ / 

(10) 

Using d^) and the fact that, by construction, the distribution function must be 
independent of Uf, one can derive a relation between the moments of the distri- 
bution function and the coefficients c^. It reads 

00 

Jdtwn{t,fi,M,z){-lnt)''~^ = c^{^,M,z); n>l. (11) 


By inverting this relation one can in principle obtain the distribution function 
from the knowledge of the set of coefficients {c^}. This observation is crucial, as 
it implies that in order to construct the distribution function it is sufficient to 
consider the perturbative series (|) in the fictitious limit /?o — *■ 00 (or n/ —00). 

To proceed, it is convenient to construct a generating function for the coeffi- 
cients c^, which we define as 

Sn{u, ^,M,z) = J2—- c^(/x, M, z) , (12) 
n=i r(ri) 

so that 

c^(/x, M,z)=i—] Sn{u, /i, M, z) _ . (13) 



du 



u=0 



This generating function is the Borel transform of the series (^ with respect to 
the inverse coupling constant 0], in the limit /3o — >■ 00. Eq. ([T^ ) can be formally 
inverted to give 



Sres{M\ z) = ^ [du Sn{% yu, M, z) exp ( - ^^"^ ) . (14) 
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In cases where the integral exists, this equation defines the Borel sum of the par- 
tial series SresiM"^ , z). In general, however, the Borel transform has singularities 
on the real w-axis, corresponding a factorial growth of the expansion coefficients 
c^. Much of the non-perturbative structure of QCD can be inferred from a 
study of the Borel transform [^-[^|. Its singularities on the negative axis arise 



from the large-momentum region in Feynman diagrams and are called ultraviolet 
renormalons. They are Borel summable and pose no problem to performing the 
Laplace integral in ([T^). Singularities on the positive axis arise from the low- 
momentum region in Feynman diagrams and are called infrared renormalons. 
Their presence leads to an ambiguity in the evaluation of the Laplace integral. In 
the following section we will discuss how this ambiguity is reflected in the integral 
over the distribution function in (|^). 

The Borel transform deflned in (|12D can be calculated by evaluating one-loop 
diagrams in which the gluon propagator in Landau gauge is replaced by 

where C is a scheme- dependent constant related to the flnite part of a renor- 
malized fermion-loop insertion on the gluon propagator. Since we have assumed 
that the expansion coefficients are dimensionless, it follows that the Borel 
transform can be factorized as 

S^{u,^^,M,z)={—^\ (16) 

where the new function S{u, z) is scale- and scheme-independent. From (0) 
and (]T^), we find that the Borel transform can be expressed in terms of the 
distribution function by the integral relation 

oo 

5r(m, /i, M,z) = J dt WK{t, /X, M, z) , (17) 



which can be inverted to give 

wn{t,^^,M,z) = ^^ J duS{u,z)l-^\ . (18) 

UQ—ioo 

The choice of Uq is arbitrary provided the integral exists. It follows that the 
distribution density wr dt depends on t, /i and M only in the combination 
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If we introduce r as a new variable and define a new scale- and scheme-independent 
function w{t, z) so that 



Q{r,z) = — I du5(u,2)r"-^ 

ZTTZ J 

UQ—ioo 

oo 



S{u,z) = JdTw{T,z)T 


oo 

C^{1,Z) = jdTw{T,z){-C -hlTY^\ (20) 


eq. ( pUD can be written in a form that makes explicit the renormalization-scheme 
invariance of the perturbative series: 

SUM\z) = Jdrw{r,z) ^^^^^ = £ [^^)\^{^., M, z) (3^^ (21) 



Note that the scheme-dependence of the constant C is such that the value of the 
coupling constant as(e'"/i^) is scheme- independent. This implies that the product 
e~'"/^AQCD is scheme-independent, where Aqcd is the scale parameter in the one- 
loop expression for the running coupling constant. We note that C = —5/3 in 
the MS scheme, C = —5/3 + 7 — ln47r in the MS scheme, and C = in the 
so-called V scheme [Q. 

As pointed out in the introduction, the integral over the distribution function 
in (|21|) is an improved one-loop approximation to the quantity S{M'^,z). It is 
instructive to compare this approximation to the BLM scale-setting prescription. 
We find 

aJre^M^) ^JaJTe^M^)\ 



SresiM, z)= dr w{t, z) ^ = 







47r \ TT 



^^Mlm) /l + A f^^^Mii^V + •••)> (22) 
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We use the symbol 



oo 



JdTid{T,z) /{t) 

(fir)) = ^ (24) 

/ dr 'u}(r, z) 



for the average of a function /(r) over the distribution w{t, z). Both the value 
of the coupling constant (ysif^BLu) ^hc parameter A are renormalization- 
scheme invariant. We observe that the first correction to the BLM scheme is 
related to the width of the distribution function. Note, however, that in some 
cases the distribution function is not of a definite sign and has no probabilistic 
interpretation. Thus, it may happen that is negative or that the normalization 
integral N vanishes, in which case (/(t)) would be ill-defined. If the distribution 
function has a definite sign, on the other hand, the quantity 

'^^^^ - I 4. J -[ln(MVA^) + (In r)r ^ ^ 

provides a measure of the rate of convergence of the perturbative series. Here 

Av = e-^/=^AQCD = e5/'3Aj^ (26) 

is a scheme-independent parameter, which coincides with the QCD scale param- 
eter in the V scheme. 

At this point it is worthwhile to point out the advantages of our resummation 
over the BLM prescription. In our approach all terms of order Pq'^u'^ in the 
perturbative series are resummed exactly, whereas the BLM scheme only resums 
the two-loop term of order j3o a"^ correctly. Moreover, our scheme is additive in 
the sense that the distribution function for the sum of two quantities is the sum 
of the individual distribution functions: wa+b = wa + wb- No such relation 
exists for the corresponding BLM scales. In particular, both A and B can have 
small BLM scales, but the BLM scale for the sum {A + B) can be large, or vice 
versa. This brings us to the most important point, which is that the size of the 
BLM scale cannot always be used as an indicator for the rate of convergence of 
a perturbative series. The size of higher-order coefficients depends on the size 
of moments of the distribution function, not only on its central value, which 
determines the BLM scale. 



3 Asymptotic behaviour, infrared renormalons 
and the separation of long-distance contribu- 
tions 

Before we calculate the distribution function for specific examples we address the 
relevance of its asymptotic behaviour for large and small values of r. Clearly, the 
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behaviour for r — oo is related to the ultraviolet properties of the perturbative 
series. In particular, only if w{t, z) vanishes faster than 1/r the integral in (^) 
converges. Otherwise, our approach provides the framework for a consistent cutoff 
regularization of the series. Performing the integral up to a value ruv corresponds 
to a hard momentum cutoff Ayy = t\j\M'^. 

More subtle is the infrared region r — > 0. As long as one stays within per- 
turbation theory one faces the problem that the integral over the distribution 
function runs over the Landau pole in the running coupling constant. One is 
forced to specify how to treat this pole, for instance by deforming the integration 
contour. In general, we may write 

/3o«s (re^M2) 1 ^ 



4:71 



^-^^ = f( ] +ri6(\nT -Ihtl), (27) 

lnMVA|. + lnr Vlnr-lnrL; ' ^ lj , k ) 



where tl = Ay/M^ is the position of the Landau pole, "P" denotes the principle 
value, and r/ is a complex parameter of order unity, which depends on the regular- 
ization prescription. This prescription dependence leads to an intrinsic ambiguity 
in the perturbative definition of S{M'^,z), reflecting the fact that in ( pTD we are 
trying to sum up a series which is not Borel summable. This is how infrared 
renormalons make their appearance in our approach. Following Refs. [p8|l-[^. 



we define the renormalon ambiguity A^rcn in the value of S{M , z) as the coeffi- 
cient of 1] and find 

— W[Tj..Z] ~ z 



. (28) 

In the last step we have used the fact that ^ 1 to expand the distribution 
function: 

w;(r, z) = woiz) t''^^ + ... for r 0. (29) 

It is thus the asymptotic behaviour of the distribution function that determines 
the size of the renormalon ambiguity. Note that A; > in order for the integral 
over the distribution function to be infrared convergent. We will see in examples 
that the power k is related to the position of the nearest infrared renormalon pole 
in the Borel transform S{u, z), which is located at mjr = k in the Borel plane. 

The appearance of infrared renormalons acts as a reminder that in (|2^ ) one is 
using perturbation theory in a regime where it is known to break down, namely 
in the infrared region. Hence, the result of any perturbative calculation in QCD 
is incomplete; it must be supplemented by non-perturbative corrections. Only 
the sum of all perturbative and non-perturbative contributions is unambiguous. 
Unlike any finite-order calculation, the representation (^) makes explicit that 
perturbative calculations contain long-distance contributions from the region of 
low momenta in Feynman diagrams. Moreover, it provides a convenient way 
to separate these long-distance contributions from the short-distance ones by 
introducing a hard factorization scale A. Thus, our approach can be used to 
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implement Wilson's construction of the operator product expansion (OPE) |]3T| in 
a literal way. Let us recall that the OPE is not designed to separate perturbative 
and non-perturbative effects, but to disentangle the physics on different length 
scales. In the calculation of short- distance corrections (the Wilson coefficient 
functions) one eliminates the contributions from small virtualities {k < A) in 
Feynman diagrams. These long-distance contributions are attributed to some 
matrix elements of higher-dimensional operators. Thus, we should write 

S,,,{M\z)= j dTw{T,z)^^^ >-+ j dTw{T,z)^^ 

A2/M2 

= S,a{M\\z) + Si^{M\\z). (30) 



Note that the factorization point r = A^/M^ corresponds to a scheme-dependent 
scale /i^ = e^\^ in the running coupling constant, with /x = A in the V scheme. 
The value of as(e'"A^) is scheme- independent, however. As long as A is cho- 
sen large enough, the short- distance contribution Ssd can be reliably calculated 
in perturbation theory and is free of renormalon ambiguities. The long-distance 
contribution must be combined with other non-perturbative corrections. Only 
the sum of all long-distance contributions is well defined. Of course, the depen- 
dence on the arbitrary scale A must cancel in the final result. This A-dependence 
can be controlled in perturbation theory by means of the renormalization-group 
equation 

A ^ 5.d(M^ A, ^) = -A ^ ^ld(M^ A, z) = -^^^^ ^ ^^(AVM^ z) . (31) 

Since we have required that the quantity S'(M^, z) be infrared safe, the long- 
distance contribution 5*^ is finite, and it is usually assumed that it is small com- 
pared to the short-distance contribution. Eq. ( pO|) allows us to quantify this state- 
ment. If the factorization scale is chosen such that A <^ M, the M-dependence of 
the long-distance contribution is again determined by the asymptotic behaviour 
of the distribution function for small values of r. We find 

.„(M^ A. .) . /V ^ . ± m) . (32) 

^ ^ 

where A(A) is of order the QCD scale, and the sign is determined by the sign 
of Wo- The long-distance contribution is indeed parametrically small, suppressed 
by inverse powers of the large mass scale M. It is, of course, no accident that 
long-distance effects appear at the same order in 1/M as infrared renormalon 
ambiguities. Both effects are exponentially small in the coupling constant and 
thus not seen in any finite-order perturbative calculation. Nevertheless, we will 
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see that in cases where there is a nearby infrared renormalon (i.e. when the power 
k is small) the long-distance contribution and the renormalon ambiguity can be 
numerically significant. One can try to estimate the size of Sid by incorporating 
certain non-perturbative effects into the integral over the distribution function, 
for instance by using a more realistic ansatz for the running coupling constant 
in the infrared region. For instance, one may assume that the coupling constant 
as(yU^) stays positive and approaches a constant for /i — > 0. It is then possible 
to perform the r-integral without encountering a Landau pole. This yields to an 
unambiguous (though mo del- dependent) result for the long-distance contribution. 



4 Heavy quark systems 

We now illustrate the formalism developed above with some quantities related 
to heavy quarks, which provide prototype examples for large-scale problems in 
QCD. The large mass scale M is provided by a heavy quark mass rriQ. Sects. ^?T| - 
|4.3| deal with the derivation of the distribution function for several quantities of 
interest. In Sect. [4.4| we present a numerical analysis of the results. 



4.1 Pole mass of a heavy quark 

We start with the relation between the pole mass mg and the (infinite) bare mass 
mo, which appears as a parameter in the QCD Lagrangian. We define 

Smiml) = ^ - 1 (33) 
^ mo 

and aim for a representation of this quantity as an integral over a distribution 
function Wm{T), as shown in (p^). In the large- /3o limit, the Borel transform 



corresponding to the perturbative series for Smimg) has the simple form 



S^{u) = 6C, (1 - u) ^^'"^^3^^^^^^^ , (34) 
where Ci? = 4/3isa colour factor. To derive the distribution function, we rewrite 

(l-«)^^^^^f^^ = ]dx]dyx-~\l-x)-^-y'~- (35) 

^ ' 

and use the second relation in (|^) to obtain 

wUr) = Qc/l^ jdyy6(T- ^^'^^'y y (36) 
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After a straightforward calculation, we find 



w,n{r) = CFr- + (l-^-]Jl + -]. (37) 



2/ V r 



The small-r behaviour of WmiT) is 



2C 

w^{t) = ^ + O(v^) , (38) 

corresponding to = 1/2 in (|29|). This behaviour is associated with the infrared 
renormalon pole at u = 1/2 in the Borel transform in (0). According to (pSD , 
the corresponding ambiguity in the perturbative series for niQ/mQ is 

,AS„U = 1^=521= = -^i^, (39) 
implying an ambiguity (AmQ)ren = (8/3/?o)Ay in the value of the pole mass 

mm- 

For large values of r the distribution function decreases as ?>CfIt-, so that 
the integral in (^) is logarithmically divergent. The divergence is removed by 
a renormalization of the bare mass. It turns out that the calculation of the 
distribution function is complicated if one chooses the MS scheme for this purpose. 
The Borel transform corresponding to the ratio vtiq jrnQimQ) is given by 



(-) = C, |6(1 - u) + e-'-" ( - ^ + R{u)) } , (40) 

where 

flM = -M"+(c(3)-^).^ + ... (41) 

is a rather complicated function. Instead, we find it instructive to consider a class 
of renormalization schemes R[r] which are more convenient for our calculation. 
We define the Borel transform corresponding to the ratio mg/mQ as 

S>) = cj6(l-„)£MI<i^-ie-| (42) 

y 1 (3 — nj u J 

and treat r as a free parameter. From a comparison of the Borel transforms 5*,^^ 
and 5*^ one finds that the relation between the so-defined mass mg and the mass 
renormalized in the MS scheme is given by 

mg _ ^ _ / _ 5\ «s ^ ^ / 2 _ lOr , 15^ fa. 



t^q{t^q) V 6/ tt 8 V 3 4 



7r 
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where = as{mQ) denotes the couphng constant in the MS scheme. It is a 
simple exercise to calculate the distribution function for the ratio rriQ/mQ. The 
result is 

w^{T) = wUr)-^&{T-en, (44) 

T 

which simply amounts to a subtraction of the high-momentum contributions, 
leaving the low-momentum region unaffected. The subtracted distribution func- 
tion falls off as for large values of r, so that the integral over r is convergent. 
With the distribution function WmiT) we compute 

A^ = i- + — ]CF = - + r, 

^2 _ 1 _ 2 2 

2.3 + | + 12C(3) + .^ (45) 

^ ^ 3 2r + l ' ^ ' 

which is all one needs to calculate the BLM scale and the parameter A defined 
in (^). Two cases are particularly interesting: 

5 

scheme Rl: r = - , 
6 




scheme R2: t = ro = \ — I ~ 2.661 . (46) 



Since the one-loop coefficient in (^) vanishes for r = 5/6, scheme Rl is similar 
to the MS scheme. The relation between the mass definitions in the two schemes 
is 

_ o ra. 



_ l + 0.208;3of—') + 0.038 /^o'f—) +... , (47) 

mQ[mQ) \ttJ \tcJ 

which is a nicely converging series. For the bottom quark the ratio equals 1.009. 
The scheme R2 is chosen such that (Inr) = 0, so that the BLM scale is given by 
/^BLM = e*"mQ. In this case one finds 

R.2 — — 2 ^3 

= 1 - 1.827^ + 0.245 /^of—) +0.006 P'J^) +..., (48) 

which has a sizeable one-loop coefficient but is again nicely converging. For the 
bottom quark the ratio equals 0.883. 

In Table |l] we give, for the case of the bottom quark, the results for the BLM 
scale in three different mass renormalization schemes (Rl, R2, MS) and for two 
renormalization schemes for the coupling constant (V and MS). We also quote the 
corresponding scheme-independent values of (^sil^^Bhydi the parameters 
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Table 1: Parameters related to the distribution function w^It) for three different 
mass renormalization schemes. We use the one-loop expression for the running 
coupling constant in the MS scheme, normalized such that as{rnl) = 0.218 for 
TUb = 4.8 GeV. 





A'-BLM 


MS 
PBLM 


O^sI/^BLm) 


A 


^BLM 


rrib / mf-^ 


0.302 rrib 


0.131 Wife 


0.55 


0.672 


0.10 


mb/mf^ 


rrib 


0.435 nib 


0.29 


4.628 


0.17 


mb/mb{mb) 


0.221 mf, 


0.096 nib 


0.73 


-4.337 


-1.19 



A and (5blm defined in (^3]) and ([25|). We note that the distribution function 
corresponding to the MS scheme must satisfy 

N = C, = \, 

^3 vr^ 

(Inr) = -— ~ -3.019, 

^ ^ 96 4 

(lnV) = ^C(3)-^ + ^^ 4.780, (49) 

as can be derived from an expansion of the Borel transform (^Oj) in powers of u. 
The large negative values of the parameters A and 5blm for the ratio nib/nib{mb) 
indicate that there are large corrections to the BLM scheme which decrease the 
value of the perturbative series. In other words, the BLM scale is too low and does 
not really represent an "average" virtuality. The reason for the strong dependence 
of the BLM scale on the subtraction scheme becomes apparent from the shape 
of the distribution functions shown in Fig. |l|. We find it most useful to show the 
product tw{t) as a function of Inr, since then the integrals (ln"r) have a direct 
graphical interpretation. The long arrows indicate the position of the BLM scale 
in the schemes Rl and R2. The small arrow shows the point r = A^/m^ for A = 1 
GeV, which will later be used to separate short- and long-distance contributions 
[cf. (pODj. In order to associate mass scales with the r- values in the figure, we 
note that ^ 

ln-^ = lnr + C, (50) 
mi 

where is the scale in the running coupling constant. In the V scheme, where 
C = 0, the point r = A^/m^ corresponds to the scale fi = X, whereas Inr = 
corresponds to /i = m^. In the MS scheme, where C = —5/3, the point fj, = nib 
corresponds to Inr = 5/3. 

We observe that the pole mass gets contributions from all momentum scales. 
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Figure 1: Distribution function for the ratio rrib/mf- in the two renor- 
mahzation schemes Rl (sohd hne) and R2 (dashed-dotted hne). The 
dotted hne shows the unsubtracted distribution function for the ratio 
of the pole mass and the bare mass. The long arrows indicate the po- 
sition of the BLM scale in the schemes Rl and R2. The small arrow 
shows the factorization point, which separates short- and long-distance 
contributions. 

and it is only the subtraction of the high-momentum tail that leads to a neg- 
ative value of (Inr) in the scheme Rl (and similar for the MS scheme). If the 
subtraction point is chosen as low as in Rl, significant cancellations take place 
between positive and negative contributions in the integral over the distribution 
function. The results is a small one-loop coefficient N in (p3|), yielding in turn 
a small BLM scale. Thus, the interpretation of the BLM scale as a "typical" 
scale in a process becomes misleading if the distribution function gets contribu- 
tions of opposite sign. Note that the situation encountered here is generic for 
quantities which require a subtraction of ultraviolet divergences. In such cases, a 
low value of the BLM scale does not necessarily imply a bad convergence of the 
perturbative series. 

The series in (|^) is interesting by itself. Since both ttiq and mQ{mQ) are 
subtracted at a large mass scale, their ratio has a well-behaved expansion in 
powers of as{frLQ) as long as r is of order unity. In fact, for the two choices 
of r discussed above the series was rapidly converging. Nevertheless, the BLM 
scale corresponding to this series exhibits a strong dependence on r. We find 
/^BLM = ^Q? with 




3 



(51) 



This function is shown in Fig. ^ The reason for the dramatic scheme dependence 
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of /iBLM is simply that the one-loop coefficient vanishes for r = 5/6. This ex- 
ample shows that a small BLM scale is not always related to a badly converging 
series. It also shows the disadvantage of non-additivity of the BLM prescription: 
As shown in Table the BLM scale for the ratio mQ/rnQ{mQ) is rather low. 
However, the same ratio can be obtained by combining the two series for mq / ttiq 
and m^/rnQ{'mQ) in the scheme R2, both of which have a much larger BLM 
scale. On the other hand, calculating mQ/mQ{mQ) by combining integrals over 
the appropriate distribution functions one obtains a unique result, in which the 
contributions from all scales are properly taken into account. 



4.2 Matching coefficients of heavy— light currents 

A convenient way to analyse hadronic matrix elements of the flavour-changing 
weak current = ^7'^(1 — 75) Q between a hadron containing the heavy quark Q 
and some light final state is to go over to an effective theory, the so-called heavy 
quark effective theory ||32|, in which such matrix elements are systematically 



expanded in powers of l/mq. When QCD is matched onto the effective theory, 
the current gets replaced by p3[ 



r ^ Ci(my q^il - 75) K + C2{ml) qv^^il + 75) K + 0(l/mQ) , (52) 

where v is the four-velocity of the hadron that contains the heavy quark, and is 
the velocity-dependent heavy quark field in the effective theory. The above form 
of the currents is correct if one uses a regularization scheme with anti-commuting 
75. The matching coefficients Ci{m^) can be calculated in perturbation theory 
by comparing quark matrix elements of the currents in QCD and in the effective 
theory. We define 

S,{ml) = Ciml) - 1 , S,{ml) = dml) . (53) 
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Table 2: Parameters related to the distribution function ^2(^)5 which is relevant 
for the ratio of heavy meson decay constants. 





A'-BLM 


A'-BLM 




xstat / -f Stat 

Jb I J b* 


0.472 mb 


0.205 mb 


0.41 2.798 7.830 0.67 



In the large-/3o limit, the Borel transforms of the perturbative series for these 
quantities are given by |^ 



S^{u) = Cf{3u^-u-3) 



2 ^r{u)T{l-2u) 



T{3-u) 



o / A r(i + u)r(i-2M) , ^, 

1 (3 — MJ 

The corresponding distribution functions can be calculated as outlined in the 
previous section. We find 



^ , , J, 7 3 ll-7r / 4 



ffi2(r) = ^|(l+T)^l + i-3-TK (55) 

The asymptotic behaviour for small values of r is 

^i(r) = -^ + 0(l), «;2(r) = ^ + 0(l), (56) 

corresponding to the infrared renormalon poles at u = 1/2 in the Borel transforms 
in (p^. If we relate the corresponding renormalon ambiguities to the ambiguity 
in the value of the pole mass [cf. (|3^) and (p9D ], we recover the relations 

(ACO„., = -H(^!!!«)=, (AC.U = |<^!!^ (57) 

12 niQ 3 niQ 

derived in Ref. The behaviour of the distribution functions for large values 



of r is ■iDi(r) ~ 1/r and W2(t) ~ 1/r^. The slow fall-off of Wi{t) leads to a 
logarithmic divergence in Ci(mQ), which must be removed by renormalization. 

As an application of these results, consider the ratio of the decay constants of 
the pseudoscalar and vector mesons B and B* in the so-called static limit, where 
terms of order Aqcd /''^b are neglected on the level of hadronic matrix elements. 



17 



In this limit, one finds 



LB 

fp^ ^ ' Ci{m 

where the eUipses represent terms not resummed in our approach. Since Ci(mQ) = 
1 + 0(1/ (3o), the distribution function is the same as for the matching coefficient 

For the distribution function W2{t) we compute 

Cf 2 



N 
(Inr) 



2 3 



2 



7 97r2 

(InV) = ^ + ^ . (59) 

The resulting value of the BLM scale and some other parameters are summarized 
in Table Since in this case the distribution function is positive definite, the 
parameter cr,- = corresponds to the width of the distribution function. The 
width is very large, about three units in Inr. This is clearly reflected in the 
shape of the function t{u2{t) shown in Fig. ^. As indicated by the position 
of the arrows, the BLM scale is rather low, close to the factorization point; the 
distribution is broad and extends well into the infrared region. As a consequence, 
the convergence of the series is bad, as reflected in the large value of the parameter 
5blm in Table ||. 



0.6 




In T 

Figure 3: Distribution function for the ratio fs^^/fB*^- The right arrow 
indicates the BLM scale, the left one the factorization point. 

In the case of the ratio mg/mg considered above, the BLM scale was low 
since the high-momentum contributions were removed, in an ad hoc way, by the 
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renormalization of the bare quark mass. On the contrary, in the present case 
there is a physical reason why the "typical" momenta in the one-loop calculation 
are much lower than the "natural" scale rrib- The flavour-changing vector and 
axial vector currents are partially conserved, implying that they do not receive 
radiative corrections from scales much above the masses of their component fields 
^ . Roughly speaking, then, the matching coefficients receive contributions from 
scales < /i < rriQ, which is in fact the behaviour reflected in Fig. ^. (Recall 
that /i = rriQ corresponds to Inr = in the V scheme, and In r = 5/3 in the MS 
scheme.) How fast the distribution falls off in the infrared region is determined 
by the location of the nearest infrared renormalon, which determines the low- 
r behaviour of the distribution function. Since in the present case the nearest 
renormalon is located at m = 1/2 in the Borel plane, i.e. at the smallest possible 
value of u, the fall-off is the slowest possible one, and thus there are substantial 
contributions from the low-momentum region. 

4.3 Matching coefficients of heavy— heavy currents 

Let us now consider a two-scale problem, namely the matching of flavour-changing 
currents containing two heavy quarks onto their counterparts in the heavy quark 
effective theory. At the so-called zero recoil point, where the two heavy quarks 
move at the same velocity, one finds 



crb^VvKrhl + 0{l/ml), 
C7^5 b ^ VA K^'l, K + 0{l/ml) . (60) 

The coefficients riv and tja take into account finite renormalization effects. The co- 
efficient r]A of the axial vector current plays a crucial role in the model-independent 
determination of | Vcb\ from B D*iiJg decays [^. 

We choose M = y/mbrric as the "natural" mass scale and define the ratio 
z = rric/mb. It is convenient to study the quantities 

Sv{mbmc,z) = r]v - I , Sv-A{mbmc, z) = r]v - r]A , (61) 
since the structure of the perturbative corrections is very similar for rjv and rjA- 



The corresponding Borel transforms in the large-/3o limit are |30, 37 



5 , , ^ r(M)r(l-2M) [2(1 + ^) 2(1 - n) - 

OvlM, z) = Gp — ^ < 1 

^ ' ' T{2-u) \ 2-u l-z l + 2u l-z 
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After a straightforward calculation, we obtain the distribution functions 

Cf< 



wv[r, z) 



l-r(l-z)2 3 (l-z^Y 

T 



2 z 
3 z 



4 



z^ 



' z 3 



TZ 



2 z 



T Z 



I 4:Z 
1 + — 



H — T z' 

4 4 T 1-z 2 1-z 

1 3 r 2 1 r 1 



wv-a{t,z) = 2Cir|l + r 

1 



4z 4^2 rl-z 2z{l-z) 

1 + Z + 

si 



1 + 



3(1-^: 

The asymptotic behaviour for small values of r is 

'i-zY 



Az / 

1 + — -Z{l + Tz)s\l 



4 

TZ 



.(63) 



2z 



(64) 



ft is associated with infrared renormalon poles at m = 1 in the Borel plane. In 
fact, a careful investigation of (|62|) shows that there are no poles at u = 1/2. The 
corresponding renormalon ambiguities in the matching coefficients are 



(Ar7y)ren 
(A?7A)rcn 



3/?o 

2 /Ay 



3/?o V"^c 



Ay' 
Ay' 



(65) 



This agrees with the results obtained in Ref. For large values of r, both 

distribution functions fall off proportional to I/t^. 
With the above distribution functions, we compute 



N 
(Inr) 

(InV) 



0.236 



3 

2' 

\ h - In^-z 

3 2 3 



In^z 



3.372, 



(66) 
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Table 3: Parameters obtained from the distribution functions for the matching 
coefficients rjv and {rjv — tja) of heavy quark currents. We use rrib = 4.8 GeV and 
= 1.44 GeV. 



















A'-BLM 




"^(/^blm) 




A 


^BLM 


Vv 


2.117 y/mbrric 




0.27 


1.059 


1.121 


0.04 


Vv - Va 


lAAby/mhrric 


O.Q28y/mbmc 


0.31 


2.069 


4.281 


0.19 



for idv(T, z), and 



2 3 

(Inr) = i + 10(2;) - 0.736, 

(InV) = —-- + Inh - (Piz) ~ 4.821 , (67) 

3 2 



for wv-AiT, z), where 



-3^1nz-6 = ^- ^ + 0{\n'z) . (68) 



We have used the value z = mc/rrib = 0.3 in the numerical analysis. The corre- 
sponding BLM scales and the values of the parameters a,-, A and (5blm are given 
in Table |^. The distribution functions are shown in Fig. ^. In both cases the 
BLM scales are comfortably large and are clearly separated from the factorization 
point, which corresponds to r = X'^/mf,mc with A = 1 GeV. We note that the 
scales /i = rrib and fi = rric correspond to Inr = ±1.204 — C, with C = in the 
V scheme and C = —5/3 in the MS scheme. The distributions fall off rapidly in 
the infrared region. Therefore, both series converge much better than in the case 
of the heavy-light current considered in the previous section. 

The physical reason for this behaviour is again related to current conservation. 
As before, the currents are conserved in the ultraviolet region, i.e. for scales 
fi ^ rrib. But in contrast to heavy-light currents, currents containing two heavy 
quarks moving at the same velocity are also conserved in the infrared region, 
i.e. for fi <^ rric- In fact, at zero recoil the anomalous dimension associated 
with such currents in the heavy quark effective theory vanishes to all orders in 



perturbation theory [l38| , p9|] . This means that the matching coefficients receive 
sizeable contributions only from scales rric < /x < nib, which is in accordance with 
the behaviour exhibited in Fig. 0. 
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Figure 4: Distribution functions for rjv (solid line) and {riv—rjA) (dashed- 
dotted line). The long arrows show the BLM scales, the short arrow 
indicates the factorization point. 



4.4 Numerical analysis 

Let us now analyse our results. For each of the quantities S{M^) considered 
in the previous sections, we compare the following approximations: the one- 
and (partial) two-loop expressions evaluated using the "natural" scale M in the 
running coupling constant, the one- loop expression evaluated using the BLM 
scale, the truncated series including the first correction to the BLM scheme given 
by the term proportional to A in (p2D, the partial resummation provided by 
the integral over the distribution function. In the latter case we use the principle 
value prescription to regularize the Landau pole in the running coupling constant. 
We define 



^1_1oop(M2) 


= N 




52-loop(M2) 


= N 




^blm(M2) 


= N 


TT 


5blm*(M^) 


= N 


as(/i|LM) 




= N 


^«,(re^M2) 





1 - (C+ (Inr 



47r 



^ //?Oas(/i|LM 



4:71 




(69) 



Note that S'2_ioop takes into account only part of the two-loop corrections, namely 
those proportional to Pq a'j.. As pointed out in the introduction, for the quantities 
considered here it is known that the remaining two-loop corrections are very 
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Table 4: Comparison of various approximations for the quantities shown in the 
first column. 





Q 

•Jl-loop 


c 

•J 2 -loop 


c 

•JBLM 


c 

•JBLM* 


c 

'-'res 


A Q 


mb 1 my- 


0.092 


0.146 


0.234 


0.259 


0.178 


0.025 


mb/mb 


U.ziy 


U.Z / i 


U.Zoo 


U.ooo 


U.oU4 


n nop; 


mb/fnb{mb) 


0.092 


0.155 


0.310 


-0.059 


0.188 


0.025 


-/"Stat / xstat 

Jb I j b* 


0.046 


0.067 


0.086 


0.144 


0.076 


0.017 


Vv 


0.020 


0.020 


0.020 


0.021 


0.023 


-0.003 


Vv - Va 


0.056 


0.065 


0.067 


0.079 


0.081 


0.007 



small. We use the one-loop expression for the running coupling constant in the 
MS scheme with A5 = 111 MeV, A4 = 150 MeV and A3 = 177 MeV, so that the 
coupling constant is continuous when one crosses the quark thresholds at = 4.8 
GeV and rric = 1.44 GeV. For reference purposes we quote that as(m^) ~ 0.218 
and as{ml) ~ 0.333. Our results are shown in Table ^. In the last column 



we give the value of the renormalon ambiguity defined in (^81) . To obtain it we 
use Po = 9 and Ay = 408 MeV, corresponding to = 3. The value of AS'ren 
should be considered as an estimate of the intrinsic ambiguity in the result for 
the resummed series S'res- 

We observe that for the ratios nib/frTbi^b) and fs^^/fB*^, where the BLM 
prescription gives very low scales and the corrections to the BLM scheme are large, 
the resummation leads to results similar to the two-loop approximation. Thus, 
higher-order corrections are smaller than indicated by the BLM prescription. 
Nevertheless, in these cases the ambiguity due to the presence of the nearby 
infrared renormalon at m = 1/2 is quite significant. For the matching constants 
r)v and r)A the resummation follows the tendency indicated by the BLM scheme, 
and the results are close to what one obtains taking into account the leading 
correction to that scheme. The renormalon ambiguities are smaller in this case, 
since the nearest infrared renormalon pole is located at u = 1. 

Another way of comparing our resummation to the BLM scheme is to define, 
for each series, a scale /x* such that the one-loop correction evaluated using that 
scale reproduces the resummed series. Hence, we write 

mb as{li.) 



/r ^ 2a^ 
fP' 3 vr 
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Table 5: Comparison of the BLM scale with the scale /i* corresponding to the 
partial resummation of the series. All values refer to the MS scheme. In the V 
scheme, the scales are larger by a factor e'l^ ~ 2.3 





rrib 


0.131 


0.195 


1.49 




nib 


0.435 


0.378 


0.87 


mb/mbimb) 


rrib 


0.096 


0.179 


1.86 


-fstat / xstat 

Jb I J b* 


nib 


0.205 


0.259 


1.26 


Vv 


y/nibnic 


0.920 


0.690 


0.75 


Vv - Va 


y/rribnic 


0.628 


0.418 


0.67 



(j){nijnib) asifi*) 
7]v — i- + - 



3 71 ' 

Vv-Va = 1^ • (70) 

In Table ^ these scales are compared to the BLM scales. Note that in the cases 
with the lowest BLM scale, namely for the ratios of masses and decay constants, 
the resummation leads to a larger scale /i* > /iblm- 

Our next goal is to obtain an estimate for the relative and absolute size of 
the short- and long-distance contributions to the various quantities considered 
above. To this end we introduce a factorization scale A = 1 GeV and evaluate 
separately the two integrals Ssd and Sid defined in (|5D|). The factorization scale 
is chosen such that the value of the coupling constant as{e'-^ X'^) / n ~ 0.25 is still 
in the perturbative regime (note that e*"/^A ~ 0.43 GeV in the MS scheme). In 
order to model the long-distance contribution, we guess a "realistic" behaviour 
of the coupling constant in the infrared region. We use a modified version of the 
running coupling constant, which exhibits freezing for /i — 0: 

We shall investigate the cases where c is adjusted so that as(0) = 1, 2 and 4 (in 
the MS scheme), and interpret the dependence of the results on c as a measure 
of the model dependence. In Table | we compare the sum S'tot = 5'sd + Sia to 
the "perturbative" resummation obtained using the one-loop running coupling 
constant regulated with a principle value prescription. We also give the results for 
the short- and long-distance contributions separately. For each case, we write the 
long-distance contribution as a power correction, | S'ld] = (A/M)'^^, and quote the 
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Table 6: Comparison of resummed "perturbative" series ^res with a model cal- 
culation of the full series including long-distance effects, S'tot = 'S'sd(A) + 5'i(i(A). 
The notation is such that the main values given for S'tot, 'S'ld and A correspond 
to as{0) = 2, whereas the corrections indicated as super- and subscripts refer to 
as{0) = 1 and 4, respectively. We use A = 1 GeV for the factorization scale. 





'-^res 


•S'tot 


5sd(A) 


^id(A) 


A(A) [MeV] 


nib / 

mb/mb{mb) 

/Stat / /Stat 
JB 1 J B* 

Vv 

Vv - Va 


0.178 
0.304 
0.188 
0.076 
0.023 
0.081 


o.243;Hi 
o.369;g:°i 

n 9^9-0.039 

U.ZOZ + Q g4g 

n 1 1Q-0.023 

090+0-001 

U.UZU_o.ooi 

0.083+o:oo5 


0.138 
0.264 
0.147 
0.057 
0.024 
0.064 


inc,-o-039 

U.lUO_,_o.046 

ine;-o-039 

U.lUO_,_o.046 
Q iQcr -0.039 

n n(S9~0-023 

U-UO^+0.029 

-0 nn4+o-ooi 

n Q1Q-0.006 

U-Uiy+o.005 


OUi_,_221 

e-r,i-186 
OUi_,_221 

dUi_,_221 

300;}^ 
154;i 
363;^^ 



value of the low-energy scale A. We observe that the long-distance contribution 
is large, as big as the short-distance one, for the ratios of the heavy quark masses 
(in the schemes Rl and MS) and for the ratio of the decay constants. The reason 
is that the nearest infrared renormalon is located at u = 1/2, leading to non- 
perturbative corrections suppressed by only one power of the large mass scale 
rrib. In general, the scales A associated with the long-distance contributions are 
typical low-energy scales of QCD. 

In Figs. ^ and ^ we show our predictions for the short-distance contributions 
to the quantities mb/rnbimb), fs/ fs*, Vv and rjA as a function of the factor- 
ization scale in the range Aq < A < M, where Aq — 0.823 GeV is the point 
where Q;s(e'"AQ) = 1. These results will become useful once non-perturbative 
calculations of matrix elements performed with a hard ultraviolet cutoff become 
available. The idea to introduce a hard factorization scale to organize the heavy 



quark expansion has been put forward recently by Bigi et al. p9[. Our approach 
provides a consistent framework to implement this proposal. Let us illustrate this 
with two examples of phenomenological importance. The first is the ratio of the 
physical meson decay constants Jb and fs* , which is related to the ratio defined 
in the static limit by 

£ xstat A 

j^ = %^ + - + Oil/ml). (72) 

jB* Jb* 

The non-perturbative parameter A can be defined in terms of hadronic matrix 
elements of dimension-four operators in the heavy quark effective theory 



which have to be estimated using non-perturbative techniques such as lattice 
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Figure 5: Short-distance contributions to m^^X) / rribirnh) (solid line) and 
{fs/ fB*){X) (dashed-dotted line) as a function of the factorization scale. 



gauge theory or QCD sum rules. These matrix elements are linearly ultraviolet 
divergent. If in the calculation of A one introduces a hard ultraviolet cutoff A in 
the same way as it was done for the perturbative calculation in the static limit, 
these matrix elements contain those long-distance contributions excluded in the 
short- distance calculation. Hence, one obtains 

p- = l + SUml A) + ^ + 0{l/ml) ^ |^(A) + ^ + 0{l/ml) . (73) 
JB* Ib* rUb 

Our resummed expression for the short- distance contribution is 



S'sd(m^,A)= J drw2(T) 



(74) 



where the distribution function W2{t) has been given in (^). The A-dependence 
in ( [73D cancels between the short- and long-distance pieces. Let us note that 
often heavy quark expansions such as ([72D are written down using dimensional 
regularization, in which case there is no clear separation between short- and long- 
distance contributions. Then the perturbative series contains an infrared renor- 
malon at u = 1/2, which is exactly compensated by an ultraviolet renormalon in 
the parameter A |^. Again, the sum of all perturbative and non-perturbative 
terms is unambiguous. 

Our second example involves the matching factor rjA for a flavour-changing 
axial vector current containing two heavy quarks. From the measurement of 
the recoil spectrum in the semi-leptonic decay B D*£ug one can extract the 
product I Kb I -^(1)5 where J^(l) denotes the value of the hadronic form factor of 
the decay at the kinematic point of zero recoil |3B]. This form factor is usually 
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Figure 6: Short-distance contributions to riv{\) (solid line) and tjaW 
(daslied-dotted line) as a function of the factorization scale. 



factorized in the form J-'{1) = r^^ (1 + ^i/ma), where the quantity Si/m^ repre- 
sents non-perturbative power corrections, which from a conceptual point of view 
cannot be distinguished from the long-distance contributions to rjA- Hence, to 
separate short- and long-distance effects properly one should again introduce a 
factorization scale, so that ?7yi(A) contains all short-distance contributions, while 
(5i/m2(A) accounts for long-distance effects. The problem is that so far all cal- 
culations of 5i/m2 were based on phenomeno logical approaches that do not ac- 
count for the A-dependence [^-[0. Therefore, we have to rely on a reason- 
able guess for the factorization scale when we combine the most recent estimate 
(^i/m^ = —(5.5 ± 2.5)% 1^ with our short-distance calculation. From Fig. |^ we 
find that 0.955 < r/A(A) < 0.975 for 0.8 GeV < A < 2 GeV, which we consider a 
conservative range of values for the factorization scale. This yields 



^(1) = 7]a{X) [1 + 5i/™2(A)] = 0.91 ± 0.03 



(75) 



which is 2% larger than the result obtained in Ref. E3], and 2% smaller than the 



value quoted in Ref. . The corresponding shift in the value of | Vch \ is at the 
level of 10-3. 



5 Correlator of light vector currents 

As an important example not related to heavy quarks, we investigate the pertur- 

bative expansion of the correlator of two vector currents in the euclidean region 

(g2 = _q2 > 

t J d'xe^'^-^ { I T{fix),fiO)} I ) = (q^q^ - q'g^ , (76) 
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where j^^ = q^y'^q. For simplicity we shall consider massless quarks. The mo- 
mentum transfer provides the large mass scale. The derivative of n(Q^) with 
respect to is ultraviolet convergent. As usual, we define the D-function 



dn(g^ 

dQ2 



(77) 



In the large-/3o limit, the Borel transform of the perturbative series for SniQ"^) is 



well known fg^, P 

Sd{u) 



?,2Cf 



E 



2 - u [A;2 
-23 



k=2 



U 



|l + - 4C(3)) w + (9 - 6C(3)) + 0( 



u 



To obtain the corresponding distribution function we start from the first relation 
in (^) and set uq = 1. This gives 



vr ^2 ^ dk 



dr e 



ir In T 



1 



1 + zr 



+ /c^ 1 — ir 



(79) 



The integral can be performed closing the integration contour at infinity; however, 
it is necessary to distinguish the cases Inr > and Inr < 0. We find 



WD[r 



7 



In r I r 

4 



8Cf 

Wo{t) = 8Ci7<jl + lnr 
where Lo(x) 



+ a + r) 



_ p 



;i + r) 

-Inr 
2 



L2(-r) +lnr ln(l + r) 



1 
r 



Inr ln(l + r-^; 



r > 1 



(80) 



Jq ln(l — y) is the dilogarithm function. The distribution 
function and its first three derivatives are continuous at r = 1, but higher deriva- 
tives are not. The asymptotic behaviour for r — > is given by 



^1) 



corresponding to the infrared renormalon pole at m = 2 in the Borel transform in 
(frSl). The location of this renormalon is consistent with the structure of the OPE 
for n(Q^), in which non-perturbative corrections appear firstQ at order 
This will be discussed in more detail below. For large values of r the distribution 
function behaves like In r/r^, so that the integral over the distribution function 
is ultraviolet convergent. 

"'^The question whether there is an infrared renormalon at u — 1 in real QCD (beyond the 
large-/3o approximation) is, however, not completely settled pq, H. 
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Table 7: Parameters obtained from the distribution function corresponding to the 
D-function. We use Ql = 2 GeV^ and Ql = (20 GeV)^. 





A'-BLM 


A'-BLM 


'^sIa'-blm) 




A 


f^BLM 


Dm 


1.628^0? 


o.rosygf 


0.40 


1.620 


2.625 


0.22 


D{Ql) 


1.628-\/gi 


0.708^01 


0.17 


1.620 


2.625 


0.03 



Using the above result, one can compute 

N = \Cf = 1, 
23 

(Inr) = 4C(3) - — ~ 0.975, 
6 

(InV) = 18 - 12C(3) ~ 3.575 . (82) 

The resulting values for the BLM scale and for the parameters cx^, A and (5blm 
are shown in Table for both a small and a large value of Q^. A graphical 
representation of the distribution function is given in Fig. |^. 




In X 

Figure 7: Distribution function for the D-function. The long arrow 
shows the BLM scale. The short arrows indicate the factorization point 
r = AVg2 for Ql = 2 GeV^ (right) and Ql = (20 GeV)^ (left). 

For completeness we also discuss the resummation for the correlator n(Q^) 
itself. To obtain it, we integrate (|77D and use the fact that SniQ"^) depends on 
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only through the running couphng constant to find 



n(g2) - n(g^; 



da. 



(83) 



where /9(q;s) = das(yU^) /d In^u is the /9-function, and Qq is some arbitrary reference 
scale, which serves to subtract the ultraviolet divergence of n((5^). To obtain the 
distribution function for n(Q^) it is sufficient to use the one-loop /3-function, 



0( \ 2 



which leads to 

in' [U{Q') - U{Ql) 



res Ql /3o 



1^7^ - — / dTWo(r) In 



(84) 



(85) 



Using an integration by parts we can bring this into the standard form of the 
distribution function representation: 



47r' 



U{Q')-U{Ql) 



In— 2 + /drwn(r) 



An 



An 



, (86) 



with the distribution function 



1 

wn{r) = — / dT'w£,(r') . 



Performing the integral gives 

wnir) = -ACf |i - Inr + (^^ - ^ Inr) t 



(87) 



+ 



L2(-t) +lnT ln(l + t) 



T<1, 



Wnir) = -4Cf |l + lnT+ (^^ + ^lnr^ i 



+ 



:i+r) 



^2(-r-')-lnTln(l + r-^) 



r > 1. 



The same result can also be obtained directly by starting from the Borel transform 
of the correlator n((5^). The asymptotic behaviour for small values of r is 



Wn{T) = -3CFT + 0{T'), 



(89) 
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Table 8: Various approximations for the perturbative contribution to the D- 
function. As previously Q? = 2 GeV^ and Ql = (20 GeV)^. 





•S*!— loop 


5'2-ioop 


5'blm 


5'blm* 


5'res 


AS'rcn 


Dm 
Dm 


0.107 
0.050 


0.125 
0.054 


0.128 
0.054 


0.156 
0.055 


0.164 
0.055 


0.006 
1.5 X 10-^ 



Table 9: Comparison of the BLM scale with the scale /i* corresponding to the full 
resummation of the series. All values refer to the MS scheme. 





A^blm/ VQ^ 




Dm 


0.708 


0.485 0.685 


Dm 


0.708 


0.621 0.878 



corresponding again to an infrared renormalon at m = 2. For large values of r 
the distribution function falls off like 1/r, in accordance with the logarithmic 
ultraviolet divergence of the unsubtracted correlator. 

Let us now turn to the numerical analysis of our results. In Table |] we show 
the various approximations to the series 5'£)(Q^), as defined in (pUj). We also 
quote values for the renormalon ambiguity 




AAcn=(A5D)ren = -r 7^ • (90) 



It is apparent that the effect of the resummation is more pronounced in the case 
where is low. In this case there are significant corrections to the BLM scheme. 
For the resummed series S'res(Q^) we define a scale [i^ so that 

Bm = ^^'^^- (91) 

TT 

This scale is compared to the BLM scale in Table ^ Finally, in Table ^ we 
evaluate the short- and long-distance contributions to the D-function introducing 
a factorization scale A = 1 GeV. 

From the asymptotic behaviour of the distribution function in (^) we con- 
clude that the long-distance contribution to the perturbative series scales like 

SM\y^) = ^,- (92) 
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Table 10: Comparison of resummed "perturbative" series S^es with a model calcu- 
lation of the full series including long-distance effects, ^tot = >S'sd(A) + 5'id(A). We 
use A = 1 GeV for the factorization scale. 





>-'tcs 


Stot 


5sd(A) 


5id(A) 


A(A) [MeV] 


Dm 
Dm 


0.164 
0.005 


U. 140+0.005 

0.005 


0.118 
0.005 


0.028;°:°°! 

(1.6 ±0.2) X 10-6 


58l;i 
713^?^ 



It should be combined with non-perturbative contributions of the same magni- 
tude. For the euclidean correlator the OPE of the current product j^{x)j'^{0) 
provides the framework for a systematic incorporation of non-perturbative ef- 
fects. At order 1/(Q^)^ these effects are parametrized by the gluon condensate 



45(] . Hence, to this order we may write 



.2 ^^ , A^(A) , 2vr 



Dm = 1 + ^sd(Q^ ^) + + ^ + 

^ 1 + 5sd(g^ A) + {a, G'){X) + . . . , (93) 

where the last equation defines the scale-dependent condensate 

{asG'){X) = {asG') + ^A\X). (94) 

Ztt 

For A = 1 GeV, we find from Table that the "perturbative" contribution to 
the gluon condensate is about 0.1 GeV^, which is of the same order of magnitude 
as the "genuine" gluon condensate {agG"^) |4^. In many practical applications 



of the OPE, and in particular in the phenomenology of the QCD sum rules, it 
is assumed that the "perturbative" contributions to the vacuum condensates are 
much smaller than the "genuine" values of the condensates and can be neglected 



[4S]. Our result (M) provides a counter-example to this assertion. 



6 Conclusions 

We have proposed an extension of the BLM scale-setting prescription, which 
resums certain vacuum polarization insertions to all orders in perturbation the- 
ory. Our approach is equivalent to performing one-loop calculations with a run- 
ning coupling constant, thereby including much of the non-trivial asymptotic 
behaviour of a perturbative series. The representation of the resummed series as 
an integral over the running coupling constant with a distribution function, as 
shown in (P^, provides an intuitive picture of the distribution of virtualities in 
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a one-loop calculation. Much insight can be gained from the knowledge of the 
distribution function. Its behaviour for large and small values of the scale param- 
eter r is related in a direct way to the ultraviolet and infrared properties of the 
series. Moments of the distribution function determine the size of higher-order 
coefficients. 

By summing an infinite set of diagrams our scheme reaches beyond pertur- 
bation theory. In particular, it provides a clear separation of short- and long- 
distance effects. In any finite-order perturbative calculation non-perturbative 
effects are implicitly present due to low-momentum contributions in Feynman 
diagrams, but are not visible as they are exponentially small in the coupling con- 
stant. Yet perturbation theory "knows" about these contributions in the form 
of infrared renormalon singularities, which make a perturbative series non Borel 
summable. This means that attempts to resum the series will lead to unavoid- 
able ambiguities. In our scheme the long-distance contributions can be explicitly 
separated, since it is possible to introduce a hard momentum cutoff in a natural 
way. The size of the long-distance contributions and their dependence on the 
large mass scale of the problem is determined by the asymptotic behaviour of the 
distribution function in the infrared region. 

Our approach offers several conceptual advantages over the ELM scheme. In 
particular, it is additive and works in cases where the BLM scale is low. We 
have emphasized that the value of the BLM scale alone cannot always be taken 
as an indicator of the size of higher-order corrections or the rate of convergence 
of a perturbative series. Except in cases where the distribution function is very 
narrow, it is better to deduce this information from the distribution function, 
which properly takes into account the contribution from all mass scales. For 
instance, the distribution function corresponding to the perturbative series for 
a quantity which requires a subtraction of ultraviolet divergences typically gets 
contributions of opposite sign, in which case cancellations may occur at one-loop 
order resulting in a low value of the BLM scale. 

The implementation of our proposal is based on techniques developed for 
the analysis of renormalon chains. The distribution function can be obtained 
from the integral relations in (^) by calculating first the Borel transform of a 
perturbative series with respect to the coupling constant in the limit of large /3o- 
We have demonstrated this for several one- and two-scale problems in QCD. We 
find that in many cases the effect of the resummation is quite significant. As 
shown in Tables |^ and ^, the difference between the resummed result and the 
two-loop approximation is comparable in magnitude to the difference between 
the two-loop and the one-loop results, in accordance with the general behaviour 
expected for asymptotic series. However, we stress that whereas the size of the 
one- and two- loop coefficients are renormalization-scheme dependent, the result 
of the resummation is scale- and scheme- independent. We have associated a scale 
fi^: with the resummed series (using a principle value prescription to regulate the 
Landau pole in the running coupling constant) and compared this scale to the 
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BLM scale. Typically, the two scales can differ by as much as 50%, which is 
significant in cases where the BLM scale is low. 

Another aspect in our analysis was to investigate the relative size of short- 
and long-distance contributions to the quantities of interest. Parametrically, long- 
distance effects are exponentially small in the coupling constant; they have the 
form of power corrections. In practice, however, they can still be sizeable in 
some cases. In heavy quark systems, for instance, non-perturbative effects are 
often only suppressed by one power of the heavy quark mass. Even for the 
bottom quark they can easily reach a level of 10%, thus being as large as one- 
loop perturbative corrections. We have introduced a hard factorization scale 
A = 1 GeV and compared the short- distance contribution to a model calculation 
of long-distance effects. We find that in some cases the long-distance effects are 
as big as the short-distance ones. Then no reliable prediction can be obtained 
based on perturbation theory alone; it is necessary to include non-perturbative 
effects. 

In Ref. [|r^, we generalize our resummation procedure to the description of 
cross sections and inclusive decay rates. In the calculation of radiative corrections 
both virtual and real gluons will have to be considered, and only the sum of their 
contributions is infrared finite p6| , Clearly, in such a situation one has to 



generalize the idea of performing a one-loop calculation with a running coupling 
constant, which was the motivation for our resummation. As a consequence, 
the "linear" form of the integral representation given in (^2]) will be replaced by 
"non-linear" representations, in which instead of the coupling constant as there 
appears a function of the coupling constant. 

While this paper was in writing, I became aware of a preprint by Beneke 
and Braun [|^, who propose the same generalization of the BLM scheme. Their 
approach is similar in spirit to the one presented here, although the formalism 
differs in technical details. I am grateful to the authors for making their results 
available to me prior to publication. 
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